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Abstract 

Quantum random walks are shown to have non-intuitive dynamics which makes them an attrac- 
tive area of study for devising quantum algorithms for long-standing open problems as well as those 
arising in the field of quantum computing. In the case of continuous-time quantum random walks, 
such peculiar dynamics can arise from simple evolution operators closely resembling the quantum 
free- wave propagator. We investigate the divergence of quantum walk dynamics from the free- wave 
evolution and show that in order for continuous-time quantum walks to display their characteristic 
propagation, the state space must be discrete. This behavior rules out many continuous quantum 
systems as possible candidates for implementing continuous-time quantum random walks. 
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I. INTRODUCTION 



Quantum random walks represent a generalized version of the well known classical random 
walk, which can be elegantly described using quantum information processing terminology 
Despite their apparent connection however, dynamics of quantum random walks are 
often non-intuitive and deviate significantly from those of their classical counterparts 
Among the differences, the faster mixing and hitting times of quantum random walks are 
particularly noteworthy, making them an attractive area of study for devising efficient quan- 
tum algorithms, including thosepertaining to connectivity and graph theory [2, 3, ij], as well 
as quantum search algorithms p, ||. 

There are two broad classes of quantum random walks, namely the discrete- and 
continuous-time quantum random walks, which have independently emerged out of the study 
of unrelated physical problems. Despite their fundamentally different quantum dynamics 
however, both families of walks share similar and characteristic propagation behavior 
Strauch's recent work [9] is the latest in a line of theoretical efforts to establishing a formal 
connection between the discrete and continuous-time quantum random walks, in a manner 
similar to their classical counterparts. 

In this paper we investigate the discretization of space in continuous-time quantum walks. 
To the best of our knowledge is the first study examining the relationship between discrete- 
and continuous- space quantum random walks. In what follows we present, in Sec. [Ill an 
introductory overview of continuous-time quantum random walks. Then in Sec. IIHI we 
provide a concise definition of discrete and continuous state space, and in Sec [IV] describe 
the way in which we model the transition of the state space from discrete to continuous. 
We present our results in Sec. [V] and show how such a transition can reduce the quantum 
walk evolution to the typical quantum wave propagation in free space. In conclusion we 
discuss the implications of our findings for the physical implementation of continuous-time 
quantum random walks. In particular we argue that existing experimental implementations 
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12l ] implicitly verify our assertion that continuous-time quantum random walks 



require discrete space. 
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II. CONTINUOUS-TIME QUANTUM RANDOM WALK OVERVIEW 



Continuous-time quantum random walks were initially proposed by Farhi and Gutmann 
j^] in 1998, out of a study of computational problems formulated in terms of decision trees. 
Suppose we are given a decision tree that has TV nodes indexed by integers i = 1, . . . , N. 
We can then define an TV x TV transition rate matrix Ti with elements 



—jij for i ^ j if node i is connected to node j 

for i 7^ j if node i is not connected to node j , (1) 
Si for i = j 

where 7^ is the probability per unit time for making a transition from node i to node j and 
for Ti to be conservative 



N 



Si = 2^1^ (2) 

Defining P as the probability distribution vector for the nodes, the transitions can be de- 
scribed by 

dPCt) 

-±> = -HPW, (3) 

for which the solution is 

P(t) = exp(-?tt)P(0), (4) 

known as the master equation. 

Farhi and Gutmann's contribution was to propose using the classically constructed tran- 
sition rate matrix Ti to evolve the continuous-time state transitions quantum mechanically. 
This involved replacing the real valued probability distribution vector P(i) with a complex 
valued wave vector ijj(t) and adding the complex notation i to the evolution exponent, i.e. 

^{t) = exp(-iWt)^(0). (5) 

Hence the probability distribution vector P(t) = |^(i)| 2 and the elements of ip(t) are the 
complex amplitudes ip(i,t) = (i\i/;(t)) where \ip(t)) is the state of the entire decision tree 
system at time t. In this quantum evolution the transition matrix Ti is required to be 
Hermitian. This formulation is not limited to discission trees and can be readily applied to 
the continuous-time quantum random walk on any general undirected graph with TV vertices. 

Figure [H shows the characteristic probability distribution \ip(t = 15) | 2 of a continuous- 
time quantum random walks on a line with TV = 160 nodes, indexed by i = —79 ... 80, and 



the initial state = 0,t = 0) = 1. For this quantum walk each node is assumed to be 
connected only to its neighboring nodes by a constant transition rate 7 = 1 resulting in a 
transition rate matrix given by 

/ -2 1 
1-210 
1-21 
1-2 
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(6) 
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For comparison we have also plotted the continuous-time classical random walk probability 
distribution P(t = 15) using the same transition rate matrix and initial condition P(i = 
0,t = 0) = 1. 



III. DISCRETE VS CONTINUOUS STATE SPACE 

In Farhi and Gutmann's treatment of the quantum walk, an arbitrary graph with TV ver- 
tices can be represented as TV position states with coordinate vectors for i — 1, 2 ... N . 
These state vectors form an orthonormal basis in the TV-dimensional Hilbert space 5, that is 
(xi\xj) = Sij and the wavefunction remains normalized. The time evolution of the quantum 
walk can be considered as continuously displacing the walker (in time) by a distance £ij from 
node i to all its neighboring nodes j at the rate 7^, where 1^ = ||^ — Xj\\ is defined as the 
transition length. Since the quantum walker can only be present at positions x\, x 2 . . . x^ ) 
we say that the state space is discrete and the walker has an infinitely narrow width. In 
other words there is no uncertainty or distribution associated with the amplitude 

i)(i,t) = (xtim) ■ (7) 

A simple example of this is the quantum walk on a line, illustrated in Fig. EK, where 
the position states are discrete nodes arranged in a line and the amplitudes i/j(i,t) are 
diagrammatically represented as narrow lines over each node. Clearly in this discrete model 
the nodes can be made arbitrarily close by making £ — ► without affecting the outcome in 
any way. 

The situation changes however when the state space is continuous, meaning that vector 
x is no longer restricted to coordinates X\ , X2 • • • Xjy and the Hilbert space is an infinite 
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dimensional continuum S' . A consequence of this is that the quantum walker's position at 
Xi can now have a finite uncertainty associated with it. This is conveniently represented as 
a distribution with a finite width Ax^ that is centered at x im Taking our previous example of 
a quantum walk on a line, this continuous model is depicted in Fig. where the amplitude 
i/j(i,t) is given by the area under the Gaussian distribution centered at x^. 

Despite the continuous nature of the state vectors in <S', we can nevertheless use this 
system to simulate the quantum walk in the discrete Hilbert space S. To implement 
this virtually- discrete state space over a continuous one, we simply ensure that the states 
1^2) . . . \xjsf) remain orthogonal by requiring that the transition length £ between all the 
nodes or position states of the walk obeys £ ^> Ax. In doing so the overlap between the 
distributions at the neighboring nodes becomes negligible, which is the case in Fig. Eb. 

What we propose here is that the relationship £ 3> Ax is in fact a necessary condition 
for the quantum random walk to retain its characteristic features. In other words quantum 
random walks require a discrete or orthonormal state space. In particular we show that for 
a continuous-time continuous-space quantum walk on a line, where the walker has a finite 
width, the evolution of the walk is reduced to the conventional quantum wave propagation, 
in the limit where the transition length approaches a continuum, i.e. £ — ► (see Fig. Eb). 

In what follows we use the finite difference approximation to construct an arbitrary 
transition rate matrix TC « — |V 2 for the quantum walk on a line. We then show that for 
£ 3> Ax the time evolution results in the characteristic quantum random walk signature 
as expected. But as the transition length £ is reduced while keeping Ax unchanged, the 
propagation behavior begins to alter until it converges to that of a quantum wave in free 
space. 

IV. MODELING THE DISCRETE TO CONTINUOUS TRANSITION 

We start with the continuous state space in Fig. [3] where the continuous position state 
vector x is formally equivalent to a variable x e K along a line. This continuous position 
space is used to construct a virtually-discrete state space quantum walk on a line by ensuring 
that the condition £ 3> Ax is satisfied. The line is then broken up to adjacent segments 
of width w ) with the position of the equidistant nodes Xi,x 2 . . .x^ at the center of each 
segment. The natural uncertainty in the quantum walker's position at Xi is given by a 
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Gaussian distribution 



Gi(x,t) = Ai{t)exp (( 



x - Xi) 2 /(2Ax) 2 ) 



(8) 



where Ai(t) is a complex phase. The amplitude of the walker to be at position Xi is then 
given by 



where L max = Xi + w/2 and L min = Xi — w/2 are the upper and lower bounds of the block 
containing the distribution. The condition w 3> Ax, guarantees that the distribution only 
finds appreciable values inside this block and is approximately zero elsewhere. 

In the way we have defined the adjacent segments, the width w is equivalent to the 
transition length £. What we want to model now is the effect of reducing £ without changing 
w ) which causes the neighboring distributions to overlap as the nodes come closer (see Fig. 
HJ). Keeping w constant is justified as the distribution within the segments is associated with 
some fundamental uncertainty in the position state which is unaffected by the change in the 
transition length. The overlap region obviously grows as the nodes get closer while the state 
space itself (i.e. the line) shrinks as a result. As I — ► and the virtually-discrete state 
space approaches a continuum, it is no longer orthogonal due to the overlapping segments, 
i.e. > and consequently amplitudes i/j(i,t) as given by Eq. [9] involve summations 

over some mutual areas. 

In order to numerically evolve the quantum walk, we represent the continuous state 
space of the walk on a finite dimensional complex vector \t of length TV x A, where A is 
the integer equivalent of £, representing the number of elements between and and 
«i, &2 • • • are elements of the vector corresponding to the node positions Xi, x 2 . . . x^ 
along the continuous line. We also introduce m as the integer equivalent of w, denoting the 
number of vector elements across each segment. 

Given an arbitrary N x N transition rate matrix Tt : in order to propagate the quantum 
walk in the quasi-continuous space represented by vector \P(t), we introduce a modified 
matrix H of size TV A x TV A. To construct this matrix, we note that a transition from the 
zth to the jth node in the discrete walk corresponds to the transition of all the elements 
within segment i to segment j as depicted in Fig. [3] and El This can be represented by a 




(9) 
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block matrix 



1-3 







\ 



B {j 



V 



(10) 



/ 





which is zero everywhere except in an m x m block 



NXxNX 



(11) 



which is centered at and and X represents the identity matrix. The modified transition 
matrix is then given by 

N N 

h = EE% ( 12 ) 

1=1 j = l 

and the continuous-time quantum random walk proceeds according to 



V(t) = exp(-iHfc)*(0). 



(13) 



V. RESULTS 



In our simulations we considered two quantum walks, both with TV = 160 discrete nodes 
(enumerated as -79, ... 0, ... 80) arranged on a line, and transition matrices Hi and H2 
which were constructed as the finite difference approximations to — |V 2 . For Hi we use the 
1 st order approximation, where 

VV(«) ~ - 1) - 2^(i) + + !))• (14) 

Here we set grid spacing /i = 1, and the time parameter in ip is implicit. Transition rates 
7ij are the coefficients of — s) divided by -2, where i — j = s and zero otherwise. 7^2 is 
constructed similarly but using the 10 th order approximation, where 

V2 ^ } ~ 2520W (W - 5) - 125^(i - 4) + (15) 

1000^(z - 3) - 6000^(i - 2) + 
42000^(z - 1) - 73766^(i) + 42000^(i + 1) + 
1000^(z + 3) - 6000^(z + 2) + 
8^(i + 5) - 125^(i + 4)). 
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For each walk we then defined the vector ^ using m — 16 which we initialized so that = 
0) is zero everywhere except for a single distribution Q { ^[x^t = 0) (Eq. [8]), corresponding 
to ip(i = 0, t = 0) = 1. The state space transition from discrete to continuous was then 
simulated by computing the quantum walk evolution for 5 different values of A = 16, 4, 3, 2 
and 1. 

In simulating the walk it is possible to construct the modified transition matrices Hi and 
H 2 as described in the previous section and compute the time evolution from Eq. [131 The 
evaluation of matrix exponentials however is computationally expensive and even for modest 
choices of TV and m, large matrices have to be stored, and the evaluation time becomes 
prohibitively long, rendering this direct method impractical. Instead in the Appendix we have 
introduce an alternative computational scheme based on Fourier analysis, which provided 
us with a dramatic improvement in the computation time. 

Figures [6] and [3 show the final probability distributions for ip(t = 15) which is obtained 
using Eq. [9j For the case A = 1 we have also plotted the analytical solution for the free 
propagation of a Gaussian wave packet [13|] 



Here we can see an almost perfect convergence of the quantum walk to the familiar propa- 
gation of the Gaussian wave packet in free space. 

Utilizing the computational scheme outlined in the Appendix, extending the above anal- 
ysis to higher dimensions becomes trivial and computationally viable. Below is an example 
where we have made the transition from discrete to continuous for a quantum walk on a two- 
dimensional mesh. The transitions follow Ti\ in the x-direction and TL2 i n the y-direction. 
The resulting probability distribution for ^2d(^ = 15) is given in Figure [8] where the quantum 
walk converges from its characteristic evolution to a symmetric Gaussian distribution. 

VI. CONCLUSION 

We have shown that the evolutionary behavior of continuous-time quantum random walks 
can be critically affected by the discretization of state space. When the quantum walker has 
a finite width, a virtually-discrete state space is constructed by keeping the position states 
of the walk well separated so that their corresponding distributions do not interfere. By 




(16) 
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allowing the position states to move closer, approaching a continuum, the walker's distribu- 
tion at neighboring position states begin to interfere which alters the evolutionary behavior. 
We have moreover shown that when TL ~ — |V 2 this can even lead to the complete collapse 
of the quantum walk behavior to that of a typical quantum wave propagation in free space. 

Beyond theoretical interest however, our analysis find important implications for the 
physical implementation of continuous-time quantum random walks. Indeed a number of 
the existing experimental schemes implicitly verify the notion that state space should be 
discrete. Du et. al [l(3] demonstrated the implementation of a quantum walk on a circle 
with four nodes, using a two-qubit NMR quantum computer. Here the four dimensional 
state space spanned by the spin states of the two qubits | ||), | jj), I IT) an d | II), is clearly 
discrete in nature and the walker has an infinitely narrow width. 

Nevertheless the assumption that the walker has a finite width is in fact a realistic one 
for many other physical systems such as a single particle in real or momentum space, where 
a localized distribution invariably arises from the fundamental uncertainty in the particle's 
position for any given coordinates. In one such proposal Cote et al [HI] described a scheme 
based on ultra cold Rydberg (highly excited) 87 Rb atoms in an optical lattice. In this scheme 
the walk is taking place in real space which is clearly continuous, but the confinement of 
atomic wavefunction to individual lattice sites amounts to a virtually-discrete state space. In 
fact Cote et al prescribe an even more conservative condition: to eliminate the atoms except 
in every fifth site (spacing 25/xm) in order to achieve a better fractional definition of the atom 
separation. Similarly Solenov and Fedichkin [j^] proposed using an a ring shaped array of 
identical tunnel-coupled quantum dots to implement the continuous-time quantum random 
walk on a circle. Given the confinement of the electron wavefunction inside the individual 
quantum dots, the authors have once again implicitly described a virtually-discrete state 
space over the continuous real space. 

APPENDIX: Computational Scheme 

Here we present a computationally efficient scheme, referred to as the Fourier-shift 
method, which not only dramatically speeds up the evaluation of matrix exponential in Eq. 
[T3l but can also be readily extended to higher dimensions when the connections between 
the nodes naturally form a multidimensional mesh. 
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Our scheme is applicable to continuous-time quantum walks on 2d regular graphs. Specif- 
ically for a graph with TV vertices or nodes, if we arrange all the nodes in a line indexed as 
1, 2, . . . , AT, we require that each node i is connected to nodes j = i — d ) i — + 1, . . . , i + 
d — 1, i + d with transition rates 7 S , where s = i — j. Applying these transitions to the zth 
node we have 

d 

H#)= £7-lW + «), (A-l) 

s=—d 

where we have made the time parameter i implicit. Using the notation 

s 

Wi = 1>(i + 8), (A-2) 

where — > represents shifting the entire vector ^ (5 nodes in the positive direction if s > 0, 
s nodes in the negative direction if s < 0, and no shift for s = 0), we can rewrite Eq. I A- II as 

d / d f A 

nm = I E 7, , (A-3) 

s=—d \s=—d / ? - 

which yields 

Hi> = 7, ^ ■ (A-4) 

s=—d 

This representation is desirable since vector shifts can be readily expressed via the Fourier 
shift theorem, which in turn greatly simplifies the evaluation of exp(— ±TLt). The discrete 
Fourier shift theorem states that 

^|T| =^Wp*(s), (A-5) 

where represents the k th element of the discrete Fourier vector with = 0, 1, 2 ... AT — 1, 
and 

p k (s) = exp(27r±ks/N), (A-6) 

where k = k for k G [0, N/2], and k = k - N for k e [N/2 + 1, AT - 1]. Then applying the 
inverse transform we have 

s 

T = ^ _1 {^W®^)}, (A-7) 
where V(s) is a vector with elements Pk{s) and ® is the direct vector product such that 
(a u 02, • • • , a n ) ® (&i, 6 2 , • • • , K) = (0161, a 2 6 2 , • • • , fl^n)- 
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We can now use the above identity to rewrite Eq. IA-41 as 

Hip = J 7 ' 1 {J 7 {ip} ® Q} , (A-8) 

where 

d 

Q=Y.^V{s). (A-9) 

s=—d 

This results in the following simplification 

U n 4> = T~ x {. ..T{T~ X {T{i)} ® Q}} <g> Q. . .} (A-10) 
= T~ x {T W <g> Q n } , (A-ll) 

where Q n denotes the direct product of n Q-vectors. Hence a polynomial expansion of Eq. 
[5] gives 



^(t) = ^C n (-iWt)>(0) (A-12) 

n=0 

= T~ x I T {^(0)} ® jr C n {-±Qt) n I (A-13) 
= T~ x {^{^(0)} ® ^(t)} , (A-14) 

where C n are the expansion coefficients and elements r^(i) of vector 1Z(t) are related to 
elements % of vector Q via the relation r^(t) = exp(— i%£). 

Computationally this representation is much more efficient than a direct evaluation of 
the N x N matrix exponential in Eq. [5l since vector Q needs to be calculated once, and for 
different values of time t we only require the evaluation of scaler exponentials exp(— ±qjt) 
followed by the action of T~ X T which can be efficiently performed using FFT. 

A numerical comparison between the direct and the Fourier-shift method shows an ex- 
cellent agreement in the resulting probability distributions, while highlighting the Fourier 
method's tremendous efficiency. Setting N = 200 and using the transition rates given by Eq. 
[T4l we were able to compute exp(— ±Hit) for t = 5 and evolve ip (initialized as a Gaussian 
distribution with Ax = 2) with a relative accuracy better than 10 -16 . Repeating this for 
N ranging between 50 and 250, we also obtained an efficiency factor td/tp as a function of 
A, where td and tp are the CPU times required for the direct and Fourier-shift methods, 
respectively. This is plotted in Fig [9] along with a quadratic fit to the data. The relatively 
large deviation of the data from the mean is mainly due to the better optimization of FFT 
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packages for arrays of certain sizes, which should be considered in the actual numerical 
implement at ion . 

We can now extend Eq. I A- 141 for application to the quasi-continuous quantum walk on 
vector ^ where each node is represented by a segment of m elements in the vector. In this 
case transition from node i to node j is equivalent to the vector elements within the ith 
segment making a transition of length A \i — j\. It is then easy to defining (TV x A)-element 
equivalents for vectors V{s) ) Q and 7£, represented by V(s), Q and 7£, and adjust the vector 
shifts by a factor A. More explicitly, the elements of V{s) are given by 

p k (s) = exp(27rite/]V), (A-15) 

where N AA, and 

Xd 

Q= E % Vis). (A-16) 

s=— Xd 

TZ(t) is then constructed as before with elements fi(t) = exp(— iup) and the time evolution 
of the quantum walk over \I/ is simply given by 

tf(t) = f- 1 {f{V(0)}®K(t)} . (A-17) 
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VII. FIGURES 
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FIG. 1: Continuous-time quantum random walk probability distribution (solid) using the transition 
rate matrix in Eq. [SJ with initial condition ip(i = 0, t = 0) = 1 and evolution time t = 15. Also 
plotted (dashed) is the time evolution of a continuous-time classical random walk using the same 
transition rate matrix, initial condition P(i = 0, t — 0) = 1 and evolution time t — 25. 
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FIG. 2: a) Continuous-time quantum random walk on a line of discrete nodes. Vertical lines 
represent the walker's probability amplitude to be at each node, d) Continuous-time quantum 
random walk on a continuous line. The nodes of the walk are a distance £ apart and localized 
distributions represent the walker's probability amplitude to be at each node with an associated 
uncertainty Ax. The condition £ ^> Ax means that the nodes are well separated and their respective 
distributions do not interfere, c) The nodes of the walk have been moved closer such that £ ~ Ax. 
Consequently the localized distributions at neighboring nodes begin to overlap and interfere with 
one another. 
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FIG. 3: A virtually-discrete quantum walk on a continuous line. The line is divided into segments 
of length w containing the localized distributions and the transition length between neighboring 
nodes is t. This continuous state space is then represented on a numerical vector \I/ with integers 
m and A corresponding to parameters w and £ respectively. A transition from node i to node j 
involves making a transition from all the elements within the zth segment to the jth segment of 
the vector. 



a lull /i; 





FIG. 4: Moving the nodes closer (i.e. t\ —> £2 — > £3 • • •) while leaving w unchanged emulates a 
move from a virtually-discrete to continuous space. This results in the overlap and interferant of 
the Gaussian distributions at the neighboring nodes. 
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FIG. 5: Approaching the nodes of the quantum walk represented on vector \I/ involves reducing A 
while leaving m unchanged which leads to an overlap region where vector elements are shared by 
neighboring distributions. A transition from node i to node j involves making a transition from 
all the elements within the zth segment to the jth segment of the vector. Vector elements in the 
overlap region are subject to the cumulative transition of the overlapping nodes. 
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FIG. 6: Convergence of the continuous-time quantum random walk probability distribution to a 
simple Gaussian for transition lengths A = 16, 4, 3, 2, and 1. The evolution is carried out using the 
coefficients of Eq. [14] as the transition rates, with the initial condition = 0, t = 0) = 1 and 
evolution time t = 15. 



17 






0.02 




0.0175 




0.015 




N 0.0125 




^ 0.01 




0.0075 




0.005 


V 


0.0025 



-75-50-25 25 50 75 
Node i 



0.02 r 

0.0175 
0.015 
0.0125 
^ 0.01 
_ 0.0075 
0.005 
0.0025 




-75-50-25 25 50 75 
Node i 



75-50-25 25 50 75 
Node i 



-75-50-25 25 50 75 
Node i 



FIG. 7: Convergence of the continuous-time quantum random walk probability distribution to a 
simple Gaussian for transition lengths A = 16, 4, 3, 2, and 1. The evolution is carried out using the 
coefficients of Eq. [15] as the transition rates, with the initial condition = 0, t = 0) = 1 and 
evolution time t = 15. 




80 -80 80 -80 80 -80 



-80 




x Node 



y Node 



-80 
-40" 

x Node 1 



80 



y Node 



80 -80 



80 -80 



FIG. 8: Convergence of a two-dimensional continuous-time quantum random walk probability 
distribution to a simple Gaussian for transition lengths A = 16,4,3,2, and 1. The evolution is 
carried out using the coefficients of Eq. [14] and [15] as transition rates in the x and y directions 
respectively, with the initial condition ^2d(^ = 0,j = 0, t = 0) = 1 and evolution time t = 15. 
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FIG. 9: The computational efficiency of the Fourier-shift method over the direct method as a 
function of the number of nodes. The solid curve represents a quadratic fit to the efficiency data 
(dotted). 
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